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Abstract. The remarkable progress in cosmic microwave background (CMB) studies over 
past decade has led to the era of precision cosmology in striking agreement with the ACDM 
model. However, the lack of power in the CMB temperature anisotropies at large angular 
scales (low-.£), as has been confirmed by the recent Planck data also (up to ^ = 40), although 
statistically not very strong (less than 3cr), is still an open problem. One can avoid to seek an 
explanation for this problem by attributing the lack of power to cosmic variance or can look 
for explanations i.e., different inflationary potentials or initial conditions for inflation to begin 
with, non-trivial topology, ISW effect etc. Features in the primordial power spectrum (PPS) 
motivated by the early universe physics has been the most common solution to address this 
problem. In the present work we also follow this approach and consider a set of PPS which 
have features and constrain the parameters of those using WMAP 9 year and Planck data 
employing Markov-Chain Monte Carlo (MCMC) analysis. The prominent feature of all the 
models of PPS that we consider is an infra-red cut off which leads to suppression of power at 
large angular scales. We consider models of PPS with maximum three extra parameters and 
use Akaike information criterion (AIC) and Bayesian information criterion (BIC) of model 
selection to compare the models. For most models, we find good constraints for the cut off 
scale kc, however, for other parameters our constraints are not that good. We find that sharp 
cut off model gives best likelihood value for the WMAP 9 year data, but is as good as power 
law model according to AIC. For the joint WMAP 9 -|- Planck data set, Starobinsky model 
is slightly preferred by AIC which is also able to produce CMB power suppression up to 
^ < 30 to some extent. However, using BIC criteria, one finds model(s) with least number 
of parameters (power law model) are always preferred. 
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1 Introduction 

Anisotropies in the temperature and polarization of cosmic microwave background (CMB) 
have information about the early universe and can be used to study the primordial fluctua¬ 
tions generated during inflation which lead to structure formation in the universe [1-6]. CMB 
anisotropies are believed to be sourced by quantum fluctuations generated during inflation 
[2, 3, 7, 8]. Thus by investigation of the CMB angular power spectrum, which completely 
characterizes CMB anisotropies, one is able to probe the physics of the very early universe. 

So far most CMB experiments indicate that CMB anisotropies are statistically isotropic 
and Gaussian and so can be completely characterized by their two-point correlations or power 
spectrum [9-11]. Although, almost all the CMB observations confirm that the six parameter 
ACDM cosmological model best fits the observed data, still there are some anomalies which 
have always been present from COBE to Planck. One of such anomalies has been the lack 
of power in the CMB-TT power spectrum at large angular scales or \ow-l [12-14]. 

Recently, [15] studied the consistency of the standard ACDM model with the Planck data 
using the Crossing statistic [16-18] . Their results indicate that the Planck data is consistent 
to the concordance ACDM only at 2-2>a confldence level and lack of power at both high and 
low €s with respect to concordance model. The low power at large angular scales can be 
attributed to cosmic variance [19], still there have been efforts to explain this anomaly by 
changing the potential of inflation held (for a comprehensive review check [20]), considering 
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different initial conditions at the beginning of inflation [21-28], ISW effect [29], spatial curva¬ 
ture [30], non-trivial topology [31], geometry [32, 33], violation of statistical anisotropies [34], 
cosmological-constant type dark energy during the inflation [35], bounce from contracting 
phase to inflation [36, 37], production of primordial micro black holes (MBH) remnants in 
the very early universe [38], hemispherical anisotropy and non-gaussianity [39, 40], string 
theory [41], loop quantum cosmology [42] etc. Although, inflationary ACDM model with 
almost scale-invariant power spectrum has emerged most successful model from the recent 
observations, it is important to note it does not uniquely confirm the generic picture of the 
universe and the generalization of primordial power spectrum having additional features like 
cut off, oscillations would be crucial in identifying specific inflationary models. 

There are two main approaches which have been followed to probe the primordial power 
spectrum of curvature fluctuations generated during inflation from the CMB anisotropies. In 
the first approach no specific model of BPS is considered and the shape of PPS is directly 
reconstructed from the data by deconvolution using different techniques [43-50]. The main 
disadvantage of this approach is that the angular power spectrum does not have all the in¬ 
formation about PPS due to nature of transformation kernel and some form of regularization 
may be needed which penalizes models which have features not desired. Maximum Entropy 
method [51] has also been applied for this purpose. Planck team has used a form of reg¬ 
ularization which penalizes any model of PPS which deviates from a straight line and has 
non-zero power at those scales at which there are no constrains [52]. The second approach 
which has been used to probe the primordial power spectrum from the CMB data has been 
to consider physically motivated models of the early universe, represented by a PPS with 
features, and estimate the parameters of those from CMB data [21, 23, 36, 53-55]. Some 
of the models of PPS are inspired from the change in physical conditions during inflation 
represented by change in the potential of inflation during slow roll or initial conditions at 
the beginning of inflation. In the present work, we consider a set of models of PPS which 
were studied in [55] and try to constrain the parameters of those from the WMAP 9 year 
and Planck data. We add a couple of other models also to our analysis. All the models we 
consider have a common feature that they all have an infrared cut off in power at large scales 
and match perfectly with the standard power law model at small scales. 

The plan of this paper is as follows: In Section 2, we give a brief outline of inflationary 
framework to introduce primordial power spectrum and its parameters and also give a short 
introduction of the models (PPS) we consider for our analysis. In section 3, we present the 
results of our analysis in the form of the best fit model parameters of PPS. We also give a 
comparison of our models using Akaike information criterion {AIC) and Bayesian information 
criterion (BIC) in this section. Discussion and conclusions of our work are given in the last 
section. 
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2 Inflationary fluctuations 


Inflation is characterized by a phase in the early universe when the energy density of the 
universe is dominated by a scalar field cj) that can be characterized by a perfect fluid with 
the diagonal components of the energy-momentum tensor given by the energy density and 
pressure density respectively [7, 8, 20]. 


and 


+ 

( 2 . 1 ) 

P4> = -y{P), 

( 2 . 2 ) 


where V{4>) is the potential energy of the scalar field. 

The dynamics of the scalar field that leads to inflation is governed by the following 
equation (in FRW case): 

'(p + mp + V'ip) = 0 (2.3) 


and 
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Slow-roll inflation is characterized by two parameters e and rj which are defined as : 
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For inflation to happen we need e << 1 and \p\ « 1. 

During inflation fluctuation 6(p in the scalar field p lead to fluctuation TZ in the spatial 
curvature which are characterized by their two point correlation function under common 
assumptions (homogeneity &: isotropy): 


(7^*(k)7^(k')) = 53(k - k')A|(A:), Vo{k) = ^ A^k), (2.8) 

where the angular brackets denote an ensemble average, 6 is the Dirac delta function 
and Vo{k) is called primordial scalar power spectrum. In the standard ACDM cosmology the 
shape of the primordial power spectrum in its simplest form can be expressed in power-law 
parameterization. This model is referred to as Power Law model and can be obtained at 
leading order slow-roll approximation of the single-inflation held [56]: 

Vo{k) = As , (2.9) 
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where n* is called spectral index (tilt parameter) and is expected to be close to 1, Ag is 
spectral amplitude, is the scalar pivot which is set equal to 0.05 Mpc“^ throughout this 
work. The scalar primordial power spectrum parameters can be calculated in terms of slow 
roll parameters (e, r/) as 




V 

~ 247r2eM|, ’ 


Ug ~ 1 + 2r] — 6e. 


( 2 . 10 ) 


In addition, to the scalar primordial power spectrum, inflation also predicts a tensor 
spectrum Vtik) due to gravity-wave (tensor) perturbations which is usually written in the 
form 

InVtik) = In At +ntln^p^ , (2.11) 

where At and nt are the tensor amplitude and spectral index respectively. In terms of slow 
roll parameter these can be written as 


At ~ 


31/ 


27r2 M|, ’ 


r = ^ ~ 16e, 
Vo 



( 2 . 12 ) 


The major contribution to Pt{k) comes form the B-mode polarization and in light of recent 
questions regarding the claims of the BICEP2 results [57-59], we will not use B-mode BICEP2 
polarization data in our analysis and therefore assume r = 0 (or Pt{k) = 0) (i.e consider scalar 
perturbations only). 

The inflation theory predict a temperature fluctuations to be statistically isotropic with 
very nearly Gaussian of zero mean, consistent with current observations. It is customary 
to represent theoretical and experiment temperature power spectrum in terms of spherical 
harmonics as [60] 

OO m=£ 

= E ^iruYUn), (2.13) 

£=1 m=—i 

where h = (0, 0) is a unit direction vector on the sky, are complex quantities and Yim{n) 
are normalized spherical harmonics. Assuming CMB fluctuations to be Gaussian distributed, 
then each aim is independent with exception equal to zero and Gaussian distributed: 


P ^mm' Ci, 


(2.14) 


where Ci is called the angular power spectrum. In practice, CMB angular power spectrum 
Ci is computed using the two-point angular correlation function 

/AT AT \ 9/-I-1 

C{hi-n2) = ( —(hi) —(h2) ) = — - CeP£{ni-h2), (2.15) 

' ' i=2 

where Pi is the Legendre polynomials. The measured angular power spectrum Ci is a robust 
cosmological probe in constraining cosmological models, the position and amplitude of the 
peaks being very sensitive to important cosmological parameters. Since Thomson scattering 
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of an anisotropic radiation field also generates linear polarization [60], there are also angu¬ 
lar power spectrum due to the polarization. The polarization anisotropies have a different 
dependence on cosmological parameters than that for temperature power spectrum and can 
provide a way to break degeneracies in various cosmological parameters. Moreover, since 
polarization data is free from ISW effect one can easily separate out the contribution of low 
CMB power due to ISW and infrared cut off. The polarization can be divided into parts that 
come from curl (B-mode) and divergence (E-mode) yielding four independent angular power 
spectra as Cj'^, Cj^, Cf^. The initial power spectrum V{k) is related to the angular 
power spectrum through 

oc j d\nkV{k)Ti^{k)Tf{k), (2.16) 

where (k) is the transfer function with X representing the CMB temperature or polariza¬ 
tion. 

One of the noteworthy outcomes from recent cosmological results, especially from WMAP 
and Planck, is the possibility of obtaining structural form of the primordial power spectrum 
[61-64], which in turn has potential to differentiate strongly between various inflationary 
models dominating early universe physics. The most commonly used primordial spectrum 
is almost scale-invariant power law during inflation which went on to produce the observed 
structure in the CMB. However, different inflationary models readily accommodate differ¬ 
ent primordial spectra with radical departures from this simple picture, especially at low 
k. Moreover, recent results of WMAP and Planck have confirmed the general picture of 
the primordial power spectrum having a suppression at low k which could not be explained 
by scale-invariant power law model. Deconvolution of CMB data strongly favors a cut off 
around horizon scale 0.00001Mpc“^ < kc < 0.0009Mpc“^ followed by a bump in a primordial 
power spectrum [62, 63]. Motivated by the fact that primordial power spectra with a cut off 
should give better likelihood than scale free power law model, we will next point out various 
primordial power spectra which have cut off at low k arising due to the physics in the initial 
phase of inflationary models. 

In the present work we consider models of primordial power spectrum (BPS) which 
suppress power at large scales (small-Zc) and agree at small scales with the standard power 
law model since our aim here is to explain the deficiency of power at large angular scales in 
CMB-TT power spectrum. By considering models with a large number of fitting parameters 
it becomes easier to fit the data, however, any method of model comparison must penalize 
models with a large number of fitting parameters. In the present work, we consider models 
with at the most three extra parameters of PPS (apart from two usual parameters Ag and 
Us)- For model selection we will use AIC and BIC'm. the next section. 
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2.1 Model 1 : Power law (PL) 


We consider the standard power law power spectrum ^‘‘Voiky characterized by two parame¬ 
ters, spectral index (rig) and amplitude Ag at some pivot scale ko. Since this model is a part 
of the standard six parameters cosmological model therefore we compare the improvement 
in the likelihood as compared to this model. Voik) is give by Eq. (2.9) and can be rewritten 
as: 


ln'P(A;) = In^ls -|- (n^ — 1) In 


(2.17) 


All the models we consider in the present work can be written as modulation over the power 
law model: 


V{k) = Vo{k)xT{k,&), 


(2.18) 


where J-{k,&) is the “modulation” part and © is a vector which characterizes the extra 
parameters. 


2.2 Model 2 : Running spectral index (RN) 

Scale dependent spectral index Ug, as characterized by an extra parameter Ug called “running 
index”, has been a part of the extension of the standard six parameter cosmological model and 
is well motivated in the inflationary framework [56, 65-67]. In the slow-roll approximation, 
ag, is second order term and is of the order of 10“^ and therefore was assumed to be zero 
for the power law model discussed in previous section, ag can be calculated in terms of slow 
roll parameters as 

~-2^-M6e7/- 24e^ (2.19) 


where the slow roll parameter ^ is related to the third derivative of the inflationary potential 
V{(p) in the following way 




— ^pi 


V'V" 


( 2 . 20 ) 


Although, larger value of ag could produce suppression of CMB power, but sizable value of 
Us will amount to violation of slow roll approximation. However, there are certain models 
[65-69] where ag can be large, while still respecting the slow-roll approximation. Therefore, 
in our analysis we also consider a model with non-zero ag. 

Current CMB observations slightly favor a non-zero running model of BPS over the 
standard power law PPS and in the present work we try to constrain the parameter ag with 
the WMAP 9 year and Planck data. We use the standard parameterization for the running 
which is given by the following equation. 


\iiV{k) 


In Ag + {ug 



( 2 . 21 ) 
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2.3 Model 3 : Sharp cut off (SC) 

This model assumes that there is sharp cut off in the primordial power spectrum at large 
scale: 


P{k) 



for k > kc 
otherwise. 


This model was considered in [48, 70] and constraints were found for the cut off scale. One 
of the interesting features of this model is that it has just one extra parameter and fits the 
data as closely as the exponential cut off model with two extra parameters which is discussed 
in Sec. (2.6). 


2.4 Model 4 : P re-inflat ionary radiation domination (PIR) 

In this model we take into account the effect of a pre-inflation radiation-dominated era which 
can lead to modulations in the primordial power spectrum [23, 24, 53]. The transition from 
a pre-inflation radiation-dominated phase to de-Sitter universe was first studied by Vilenkin 
& Ford [53]: 

V{k) = ^ ^ |e-2*^(l + 2iy) -l-2y^\\ (2.22) 

4y4 

where y = k/kc- The cut off scale kc is set by the Hubble parameter and is proportional 
to k‘^. Here, current horizon crosses the horizon around the onset of inflation. This model 
produces cut off followed by a bump like feature in the primordial power spectrum. 


2.5 Model 5 : Pre-inflationary kinetic domination (PIK) 


We also consider a model given in [21] which also produces cut off due to possible existence 
of a kinetic stage in the pre-inflationary era, where the velocity of the scalar field was not 
negligibly small. In order to affect the low-.^ multipoles, this stage should occur very close 
to the beginning of the last 65 e-fold period of inflation. If scales corresponding to the 
current horizon have exited the horizon around the onset of inflation then this could cause 
a significant drop on the large angular scales of the primordial spectrum. Here the inflation 
potential is quadratic 

V(</-) = (2.23) 

with initial conditions given by = 18Mp, {d(l)/dt)in — —m(p4>in- The form of primordial 
perturbations for pre-inflationary kinetic domination model can be expressed as 

V{k) = ^k\A-B\^, (2.24) 

where 


g ik/ 




f Hjnf 

V k 
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B = 


[■^(2) ^ k ^ 

° V2i^*n/y 


f Hjnf 

V k 


n 


( 2 ) 



(2) (2) 

i/j„/ denotes the Hubble parameter in physical units during inflation, T-Lq ' and T-Li ' denote 
the Hankel function of the second kind with order 0 and 1, respectively. Here the cut off is 
proportional to k^. However, this model has a scale invariant primordial power spectrum for 
large k and therefore is strongly disfavored by the current data, despite producing low CMB 
power. 

If we consider that quantum fluctuations originate in the Bunch-Davies vacuum as is 
considered in [23] also then the primordial power spectrum can be rewritten as: 


( 1 \ ns —1 it 2 

_) (2.25) 

with 

As = A',^ko\Aiko)-B{ko)\^. (2.26) 

where ko is the pivot scale. 

This model also has one extra parameter Hinf which we constrain from WMAP 9 + 
Planck data. 


2.6 Model 6 : Exponential cnt off (EC) 

The primordial power spectrum which has lesser power at low k can also be approximated 
by imposing an exponential cut off at /c < kc [21, 70-72] on the power law model Vo{k)'. 


V{k) 


Vo{k) 



(fc/fcc)“ 


(2.27) 


where a is a measure of the steepness of the cut off. On small angular scales, this param¬ 
eterization behaves as simple power law model and qualitative features in the CMB power 
spectrum that determine the constraints on the cosmological parameters are not affected 
except at the \ow-i multipoles. 


2.7 Model 7 : Starobinsky (SB) 

Another model which predicts a step-like feature in V{k) was proposed by Starobinsky [54], 
which assumes that there is a sharp change in slope of potential of the scalar field V^cj)) at 
certain (j)Q which controls the inflationary stage. The general form of scalar field potential 
which has a rapid change for such cases can be expressed as 


v{4>) 


Ho + A+ (0 - (/>o) 0 > 00 
Vo -k A_ (0 - 0o), 0 < 00 ’ 


(2.28) 


1 


Note that the validity of imposing Bunch-Davies vacuum in Kinetic Domination regime at all scales is 


questionable but this assumption does provide a specihc form of model PPS to study. 
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where A- and are model parameters assumed to be greater than 0. It can be found that 
if the width A(j) ~ (cj) — (j)Q) of the singularity is small enough then the resulting adiabatic 
primordial spectrum is non-flat around the point kc which can be expressed analytically in 
terms of transfer function applied on any underlying power spectrum: 


V{k) = ro{k)V\y,A), 


(2.29) 


where the transfer function is given by [73, 74], 


V\y,A) = 


1 + 


9A2 




+ 


3A 


A + 3A- 


3A\ 1 


y j y 


■cos (2?/) 


+3A ( 1 - (1 3A)^ - - sin(2?/) 

y^ y^ J y 


(2.30) 


Here y = k/kc and A = —■ In this model we have applied transfer function on simple 
power law model Vo{k). kc determines the location of the step and has no effect on the shape 
of the spectrum besides the overall normalization. For R = A^jA- < 1, there is a sharp 
decrease of spectrum followed by a bump at small k with large oscillations and a flat upper 
plateau on small scales (see also [75] ). For R> 1 this picture is inverted and has a step-down 
like feature. Here, in this model for large k the modulation term becomes close to 1 and we 
get simple power law model. 


2.8 Model 8 : Starobinsky cut off (SBC) 

As discussed in the previous section, Starobinsky’s transfer function can be imposed on 
any class of primordial power spectrum. Here, in this model, we superimpose Starobinsky 
modulation on an exponential cut off spectrum (model 6) with a adjustable sharpness of the 
cut off: 


V{k) 


Vo{k) 


1 _ 


V^y^A), 


(2.31) 


where T>‘^{y,A) is the transfer function of the Starobinsky feature described in the previous 
section and e sets the ratio of the two cut off scales involved. This model has both exponential 
(sharp) cut off and a Starobinsky model bump like feature in the power spectrum. Previously 
Sinha & Souradeep [55] have found that this model provides best likelihood value among the 
wide range of models discussed here. For simplicity, we reduce the degree of freedom of this 
model by hxing e=l. We found that this parameterization does not affect final results. 


3 Methodology and parameter estimation 

We employ Monte Carlo Markov Chain (MCMC) analysis to estimate the parameters of PPS 
models we consider for our study and use publicly available code COSMDMC [76, 77] for this 
purpose. COSMQMC uses publicly available code CAME [78, 79] for computing angular power 
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spectra of CMB anisotropies following a line of sight approach which was given in [80]. 
COSMDMC uses the likelihood code provided by the WMAP and Planck team for computing 
the likelihood. 

WMAP 9 parameter estimation methodology is given in [10] which is not very different 
than what was outlined in [81]. WMAP 9 year likelihood code does not need any extra 
parameter and computes the likelihood at low and high .^’s differently for the temperature 
and polarization data. At \ow-^ [I < 32) TT likelihood is computed from the angular power 
spectrum estimated using Gibbs sampling and at high-.£ {I > 32) TT likelihood is calculated 
from the angular power spectrum estimated from an optimum quadratic estimator. For 
polarization, high-£ [I > 23) TE, EE and BB likelihoods are computed using MASTER and 
low-£ (/ < 23), TE, EE and BB likelihoods are computed in the pixel space. 

Apart from WMAP 9 year temperature and polarization data, we also consider Planck 
temperature data for our analysis which is also publicly available. Planck likelihood code 
(discussed in [14] and downloadable from [82]) has different modules to compute likelihood 
at low and high £’s. Planck likelihood code also computes likelihood for low-^ polarization 
data which it uses from WMAP 9 year data, however, we do not use that. We consider only 
modules which compute TT-likelihood at low and high At high-£ (up to £=2500), Planck 
likelihood code uses a code named CamSpec which has 14 extra parameter to take care of 
foreground and other systematic. At low-£ (£ < 49) Planck likelihood code nses COMMANDER. 

Since we perform joint WMAP 9 + Planck analysis of CMB power spectrum, it is im¬ 
portant to discuss about consistency between WMAP 9 and Planck data. Within in the 
context of the ACDM model, it has been found that the values of some cosmological param¬ 
eters like Hq obtained from Planck measurements are significantly different from WMAP 9 
measurements. It was shown in [83] that WMAP 9 angular power spectrum is about 2.6% 
higher at very high significance level at low £’s, however, no significant bias was fonnd at 
high £’s. Similarly, Planck team has reported the the inconsistency of the angular power 
spectrum at multipoles £ < 40 [14] and 2% difference of angular power spectra near the first 
first acoustic peak [84]. It was also shown by [85] that although best fit model to Planck data 
was consistent with WMAP 9 year data, but WMAP 9 best fit was found to be inconsistent 
with Planck data at 3(t. Recently [86] revisited the analysis of the WMAP 9 data and they 
also found fonnd 2.5% difference in WMAP 9 and Planck spectra at £ > 100 at 3-5a level. 
Although, some level of the inconsistency between cosmological measurements was found 
arising from the Planck 217 x 217GHz detector [87], but there still remains significant ten¬ 
sion. The tension between WMAP 9 and Planck data could be due the different systematics 
present in the WMAP 9 and Planck data or it could signal the the failure of AGDM model 
which could have far reaching consequences. 

We modify GAME and COSMDMC so that the extra parameters of the PPS models can be 
incorporated and use priors as are given in Tab. (1). Since for running COSMDMC we need 
covariance matrices also apart from prior range therefore we generate covariance matrices 
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Parameter Name 

Symbol 

Prior Ranges 

Baryon Density 


0.005-0.1 

Cold Dark Matter Density 


0.001-0.99 

Angular size of Acoustic Horizon 

e 

0.5-10.0 

Optical Depth 

T 

0.01-0.8 

Scalar Spectral Index 

Us 

0.5-1.5 

Scalar Amplitude 

loglO^^As 

2.7-4.0 

Hubble Parameter at Inhation 

H,nf (Mpc-1) 

lo-^-io-"^ 

Running Index 

(Xg 

-1-1 

Cut off Parameter 

kc (Mpc“^) 

0.0-0.01 

Cut off Steepness Parameter 

a 

1.0-15.0 

Starobinsky Parameter 

A 

0.0-1 


Table 1: Uniform prior used in parameter estimation. 


from a few trial runs. 

The cosmological parameterization has been carried out by using the six basic param¬ 
eters (baryon density cold dark matter density Thomson scattering optical 

depth due to reionization “r”, angular size of horizon “0”, spectral index “n^” and scalar 
amplitude “InlO^^A^”) along with the parameters which describe the features in the PPS 
i.e., kc, a, A etc. 

Apart from the standard six cosmological parameters, we keep the values of the rest of 
the cosmological parameters constant. We have hxed the sum of physical masses of standard 
neutrinos ‘V”=0.6 eV, effective number of neutrinos “Ag//” =3.046, Helium mass fraction 
“YHe”=0.24 and the width of reionization 0.5. For the case of WMAP 9 year + Planck data 
all the nuisance parameters of CamSpec where hxed to the standard values given in [11, 14]. 

We perform a Markov Chain Monte Carlo analysis to determine the values of the model 
parameters that provide the best ht to the observed data from Planck and WMAP for CMB 
power spectrum. Getdist was used with the chains generated by CDSMDMC to produce 2D 
contours and plots of the marginal posteriors. 

3.1 Best fit parameters 

We present the results of our analysis in terms of the best ht values and their mean values 
with l-cr errors (when possible) for the parameters of the PPS models which characterize 
the primordial power spectrum. Since we hnd that the values of the rest of the cosmological 
parameters are within acceptable range and do not show any interesting correlation with our 
model parameters of PPS, so we do not present estimates for those here. 

We present the estimates of model parameters for the WMAP 9 year and WMAP 9 year 
+ Planck data in Tab. 2 with the values of -21og likelihood (or From the table we can 
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see that all the models we consider give better fit to the data than the standard power law 
model. However, the improvement is marginal. We also present the ranking of the models 
later in this section. 



WMAP 9 

WMAP 9-hPlanck 

Model 

Parameter 

Best Fit 

68% Limit 

= -2 log/: 

Best Fit 

68 % Limit 

X^ = -2 log/: 

1 (PL) 


7558.0160 

15382.9400 

2 (RN) 

a 

-0.012 

-0.012± 0.022 

7557.7340 

-0.009 

-0.009 ± 0.006 

15380.6580 

3 (SC) 


2.9149 

2.3597± 0.9809 

7555.6080 

3.0449 

2.5653±0.8250 

15378.2840 

4 (PIR) 


0.3910 

0.5909±0.5324 

7557.9100 

0.3941 

<0.4296 

15382.1560 

5 (PIK) 


2.0846 

2.07836 ± 1.0052 

7556.1900 

2.1485 

2.0934± 0.8973 

15380.0950 

6 (EC) 


2.9244 

2.4876±1.1406 

7555.6700 

2.9780 

2.7752± 0.9237 

15378.6420 


a 

7.6167 

8.1354[NL] 


9.22328 

8.2764[NL] 


7 (SB) 

10‘^kc 

1.4724 

< 12.83404 

7556.1760 

14.641 

< 14.5739 

15375.7760 


A 

0.3893 

< 0.2583 


0.0558 

0.0696±0.0667 


8 (SBC) 

lO'^fcc 

3.1313 

< 0.1839 

7555.7640 

2.9149 

< 0.23354 

15378.3460 


a 

12.502 

8.022 [NL] 


12.6627 

8.1238[NL] 



A 

0.0037 

<0.4509 


0.055492 

<0.2896 



Table 2: The best fit and mean values of the extra parameters of PPS models we consider for the 
WMAP 9 year and joint WMAP 9 year and Planck data. We hnd good constrains on the cut off 
scales kc, however, our constrains on other parameters are poor. We were able to put upper limit on 
the Starobinsky parameter A but no limit (NL) was found for the exponential cut off parameter a. 


One of the features which all of our models (apart from PL) have common is a cut 
off at large scale we have and characterized by a scale kc- We find that cut off around 
1.40 X 10“'^ — 3.15 X 10“^Mpc“^ for most of the models discussed here, except for model 

(4) which predicts much smaller value of kc for both data sets and model (7) which predicts 
much larger values of kc for joint WMAP 9 year + Planck data set. Fig’s. (1), (2), (3), (4), 

(5) , (6) and (7) show the marginal one-dimensional posterior distributions and 2D contours 
at 68% and 95% CL of the parameters describing PPS models discussed in the work. For 
models (3), (6), (7), (8), we were able to obtain good bounds on the kc for the WMAP 9 + 
Planck data as is clear from their two dimensional joint probability distributions. Fig. (8) 
shows the WMAP 9 -|- Planck best fit primordial spectra for the range of models discussed 
in this work. Note that all the models we consider here have cut off at large scale k < kc 
and matched with the standard power law model at large k. Fig’s. (9) &: (10) show the 
corresponding angular power spectra obtained using best fit values of PPS parameters 
and other standard cosmological parameters. 

3.2 Model comparison 

There are many methods for model comparison and have their own advantages and disad¬ 
vantages. Two of the most common methods are Akaike Information Criteria or AIC [88] 
and Bayesian Information Criteria or BIC [89]. BIC is considered more powerful than AIC 
since it explicitly takes into account the number of data points when penalizing the models 


- 12 

















Figure 1: Two dimensional joint probability distributions and one dimensional marginal probability 
distribution for the parameter 10®As,ns and as of the primordial power for RN model for WMAP 9 
+ Planck data. 


with more parameters. AIC has been most commonly used [74, 90] and is defined in the 
following way: 


AIC = -2\iiC{d\e) + 2k, (3.1) 

where C{d\6) is likelihood, d data vector, 0 G is the parameter vector and k is the number 
of parameters. The best model is one which have minimum value of AIC. 5/C is also defined 
in the similar way: 

BIC = -2 In /:(d|0) + A: In iV, (3.2) 

where N is the number of data points. 

The second terms in the RHS of equation (3.1) and (3.2) work as regularization functions 
and penalize models with more number of parameters like any other regularization function 
does and can be related to Maximum Entropy Method [51]. BIC is more conservative than 
AIC since it puts stringent penalty for models with more fitting parameters for InA^ > 
2 i.e., the penalty term exceeds 2k. It is well known that AIC minimizes the Kullback- 
Leibler divergence between the estimate from a candidate model and the true distribution 
and BIC selects a model that maximizes the posterior probability distribution. For AIC, 
AAIC = AlCi — AlCmin represents preference of model i over the the best fit model. Models 
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Figure 2: Same as in Fig. (1) for the model 3 (SC). From the figure it is clear that the cut off scale 
is quite well constrained from the data. 


with AAIC < 2 have substantial support, models with 4 < AAIC < 7 have considerably 
less support and those with AAIC > 10 have essentially no support compared to best fit 
model [91]. ABIC = BICi — BICmin represents the preference of best fit model ( i.e 
model with minimum value of BIC (BICmin)) over model i. ABIC values of ABIC < 2, 
2 < ABIC < 6, 6 < ABIC < 10 and ABIC > 10 represent weak, positive, strong and very 
strong support for best fit model respectively [92]. However, in order to calculate BIC, it 
is not straightforward to find out the number of independent data points (N) for WMAP 
9 year + Planck data as it is a correlated data set. But it is clear that as the number 
of data points is quite large, BIC becomes highly sensitive towards the number of extra 
parameters. For a sensible estimates of N, we can assume 1168 (WMAP 9 Master TT £’s) + 
777 (MWAP 9 Master TE ^’s) + 1170 (WMAP 9 TT/TE/EE/BB low-^ chi2 pixels) + 2499 
(Planck CamSpec+commander £’s) = 5614 independent data points for WMAP 9 + Planck 
data set. Therefore, we get for WMAP 9 data In(A') = 8.04 and In(A') = 8.63 for WMAP 9 
+ Planck data set. 

We present AAIC and ABIC values for the models we consider in Tab. (3). We find 
that for the WMAP 9 year data, sharp cut off model (SC) gives the best (minimum) value 
of AIC but is as good as power law model, however, for WMAP 9 year + Planck data 
Starobinsky model (SB) has the lowest value of AIC which we believe is this due to the 
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Figure 3: Marginal one dimensional probability distributions and two dimensional joint probability 
distribution for the model 4 (PIR) which also has just one extra parameter for the WMAP 9 + Planck 
data. For this model the constraints on kc are not that good as we have for model 3 (SC) which is 
expected since the role of kc in this case is slightly different. 



WMAP 9 

WMAP 9+Planck 

Model 

AAIC 

ABIC 

AAIC 

ABIC 

1 (PL) 

0.408 

-5.632 

3.164 

-10.036 

2 (RN) 

2.126 

2.126 

2.882 

-3.718 

3 (SC) 

0.000 

0.000 

0.508 

-6.092 

4 (PIR) 

2.302 

2.302 

4.380 

-2.22 

5 (PIK) 

0.582 

0.582 

4.319 

-4.281 

6 (EC) 

2.062 

8.102 

2.866 

2.866 

7 (SB) 

2.56 

8.608 

0.000 

0.000 

8 (SBC) 

4.156 

16.236 

4.570 

11.17 


Table 3: This tables shows AAIC and ABIC for the different models for WMAP 9 and WMAP 
9+Planck data. For WMAP 9 data we find that the sharp cut (SC) model (model 3) gives lowest AIC, 
however, for WMAP 9 year + Planck, Starobinsky model (SB) model (model 7) gives the lowest AIC. 
The reason behind model 7 being preferred by WMAP 9 year + Planck data is that it suppresses 
power at higher angular scales i < 30. However, BIC always prefers power law model over the cut 
off models because of least number of parameters. 
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Figure 4: One dimensional marginalized probability distribution and contour plots for the parameters 
of PPS for model 5 (PIK) with WMAP 9 + Planck data. 


fact the this model gives suppression of power up a higher values of I as required by the 
Planck data. Using BIC rule for model selection, one finds that power law model is always 
favored over the cut off models. Moreover, among cut off models, model with lesser number 
of parameters is always favored. It is also worth mentioning that model (8) has a larger value 
of despite having more parameters and therefore is disfavored by the current data. 

4 Discussion and conclusions 

The observed CMB power spectrum is in striking agreement with the standard ACDM model 
with almost scale-invariant adiabatic fluctuations produced during the inflationary epoch. 
However in such studies, some anomalies have been observed such as low CMB power on 
large angular scales. It has been observed that inflationary epoch cannot be well described 
by simple form of scalar power spectrum based on the smooth slow roll approximation and 
the presence of the cut off in the primordial power spectrum is essential for extension of this 
simplistic picture in order to explain low CMB anomaly. In this work, we have explored 
different parameterizations of inflationary driven primordial spectra which have cut off at 
large angular scales so as to describe low CMB anomaly. 

We have analyzed the complete CMB data sets of WMAP 9 year and Planck. We 
perform a Markov Chain Monte Carlo analysis to determine parameters that provide the best 
fit to the data for the CMB angular power spectrum. We find that primordial power spectrum 
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Figure 5: We find that the likelihood for WMAP 9 + Planck data is not very sensitive for parameter 
a for model 6 (EC) so we have poor constrains (any value of a > 5 value is as good as any other 
value). However, for this model we also obtain good constrains on the cut off scale kc- 


with cut off, in general, leads to improvement of the likelihood and marginal preference for 
a non vanishing cut off scale of kc in some cases. Due to the large variance in the CMB 
temperature at low multipoles, we could only place weak constraints on some parameters of 
our model like a, however, our constrains on the cut off scale are fairly good. 

In order to quantify the significance of the fits we have used Akaike information criterion 
and Bayesian information criterion. We find that for the WMAP 9 data, among various 
models discussed here, model (3) which has sharp cut off provides best likelihood value, but 
is as good as power law model as per AIC. For the WMAP 9 year + Planck data set, we 
find only Starobinsky’s model (7) is able to explain suppression up to multipoles £ < 30 as is 
indicated by the much large value of the kc = 14.64 x 10“'^ Mpc“^ and is slightly preferred 
over other models as per AIC. Although, Starobinskys model improves the fit in the Planck 
CMB power spectrum in the region .^ < 30, the produced suppression is still not enough and 
there is some scope of improvement in the fit. However, using BIC, one finds that power 
law model is always preferred over all cut off models. 

It is also important to note that the power suppression in the CMB anisotropy is cur¬ 
rently a subject of intense debate. The CMB suppression could be caused by other mecha¬ 
nisms and therefore, it is important to improve and develop current techniques which allow 
us better understanding of CMB suppression. There are also models such as [93-95] which 
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Figure 6: For model 7 (SB), we find upper limits for the fitting parameters kc and A for the WMAP 
9 + Planck data. 

have oscillations superimposed on primordial power spectrum which also provide better fit to 
CMB power spectrum, but in these models the cut off is not evident as is preferred by current 
data and oscillations cover over whole range of k (unlike the models discussed here, where 
oscillations die out after the cut off). Finally, we conclude that the present motivation for the 
low CMB power at small I with an infrared cut off is very high and there is significant scope 
in improving the estimates of power suppression on the basis of modeling of the primordial 
power spectrum with the upcoming and future data. 
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Figure 7: Like model 6 (EC), we find for model 8 (SBC) poor constraints for a, however, we find 
good constrains for the parameter kc and A for the WMAP 9 + Planck data. 
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Figure 8: The best fit primordial spectra for the models consider for our analysis using WMAP 9 + 
Planck data. Note that all the models we considered in this work have cut off at large scale k < kc 
and matched with the standard power law model. All the models we consider give better likelihood 
than the pure power law model. 
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Figure 9: The best fit angular power spectra for the models of PPS we consider using WMAP 
9 + Planck data. The observed data points for WMAP 9 + Planck data are also shown by red dots 
and blue triangles respectively with error bars. 



Figure 10: Same as in Fig. (9) at low-£. This figure shows that the model 7 (SB) gives power 
suppression up to a large values of I so fit the Planck data better than other models. 
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